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ABSTRACT : 

An error in Westergaard' s equation for a certain class 

* 

of plane crack problems, originally pointed out by Sih, is 
briefly discussed anew. The source of the difficulty is 
traced to an oversight in an earlier work by MacGregor, upon 
whose work Westergaard based his equations. Several examples 
of interest, illustrating the consequences of the necessary 
correction to these equations are given. 
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INTRODUCTION 


The Westergaard equations, which apply for a certain class 
of plane problems in linear elasticity, were shown to be gen- 
erally incorrect by Sih in 1966, [1]. Specifically, by use of 
the well known Goursat-Kolosov complex representation of the 
plane problem, it was shown that the stress and displacement 
field equations appropriate to the restricted class of problems 
alluded to above include a real constant term which is lacking 
in the Westergaard equations. 

In this paper the constant term which, according to Sih's 
analysis, should be appended to Westergaard ' s equations, is shown 
to be the result of an oversight in a lesser known work of 
MacGregor [2] , upon whose work Westergaard based his formula- 
tions [3] . The consequences of the corrected equations are 
then demonstrated for several familiar plane crack problems, 
and for the approximate plane crack-tip stress and displacement 
field equations . 

The problem of the centrally cracked strip of finite width 
loaded unaxially in uniform tension is also discussed. A 
Westergaard type stress function is introduced which provides an 
approximate closed form solution. This approximate solution has 
the merits of yielding the Fedderson secant formula for the crack- 
tip stress intensity factor, and for providing an analytical 
expression for the crack opening displacement which closely 
matches experimental data and which is a considerable improve- 
ment over the calculation first introduced by Irwin [4] . 



MODIFIED WESTERGAARD EQUATIONS 

In MacGregor's complex characterization of the plane 
problem (omitting body force) the holomorphic functions 

J(z) = 0(x,y) + i fi(x,y) 

H(z) = 0 Q ( x , y ) + i fl Q (x/y) (1) 

are introduced together with their derivatives 

i J'(z) E i W(z) = $(x,y) + i V(x.,y) 

H'(z) E - K (z) = - V(x,y) - i II(x,y). (2) 

The bi-harmonic Airy stress function U(x,y) is represented 

as a linear combination of the single-valued harmonic 

functions 8 and 0 by 

o 

U(x,y > = y 0 + 0 Q = U(z,z) = Im[iJ {zf) ] 

+ Re[H(z)] (3) 

The complex representation of the plane stress field is 
then readily shown to be 

a xx = 2$ + y + = 2 Re [ iW ( z ) ] - ylm[iw'(z)] + Re[K'(z)] 

°yy = ~ Y - |- = + ylm[iw" (z) ] - Re[K"(z) ] 
a xy = ~' ? " y fx + fy = ~ i^tiWCz)] - yRe[iW"(z)] 

- Im[K"(z)] (4) 
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For the restricted class of plane problems for which 
a X y = 0 a t all points along the line y = 0, which includes 
plane crack problems in which the internal crack (or cracks) 
is situated along the x axis and where the applied loads 
are symmetrically ' located with respect to the crack plane, 
it follows from (4) that 

¥ - = Im[iW (z) + K"(z)] = C. (5) 


Consequent ly 

3x 9y 3x 

from which it necessarily follows that 


i* „ __9 ,3V, 

3x “ 3x l 3x ; 

/ 

which is satisfied in the most general sense if one chooses 


or 


Re [ iW ( z ) + K' {z) ] = - A ( 7 ) 

everywhere. Here A is a real constant. The oversight in 
MacGregor's work rests in the fact that A was omitted or, put 
another way, was necessarily presumed to be zero. Substituting 
equations (5) thru (7) into (4) and introducing 

Z(z) = i W(z) . . (3) 
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one obtains > 

G xx = Re [ Z ( z ) ] - y Im [ Z *" ( z ) ] - A 
a yy = Re; [ Z ( z)] + y I mfz"(z ) ] + A 

° xy = yRe[z'(z) ] . ( 9 ) 

which are the equations obtained by Sih when Z(z) 5 2<p'(z). 


Because the stress components are required to satisfy given 
boundary conditions the constant A will in general depend or 
the manner of the applied loading and will vanish only for 
rather special loading conditions . 

The displacement field equations must likewise be cor- 
rected. In the Goursat-Kolosov representation the displace- 
ment field is specified by the well known form [S] 

/ 

2y(u+iv) = kcJ) ( z ) - z <p > (z) - \pJzY . (10) 

where u(x,y) and v(x,y) are respectively the x and y com- 
ponents of the displacement vector, y = E/2(l+v) is the shear 
modulus, E and v are Young's Modulus and Poisson's Ratio 
respectively, and k = [3-v/l+v] for plane stress and k = [3-4v] 
for plane strain. The holomorphic functions <J>(z) and \p(z) can 
be shown to be related to those introduced in equations (1) , 

(2) and (3) by the relations 

iW (z) = Z (z) = 2(f>' (z) 

H(z) = X(z) + zcf(z) 


4 



( 11 ) 


H'(z) - - K ( z) = p(z) + <p(z) + zip' (z) 

X'(z) ~ ip( z) . . . . 

Adding to (10) its complex conjugate in the case of plane 
stress, one obtains 

Eu = (3-v) Re [<f> (zj ] - (1+v) {xRe'[<j) " (z) ] + ylm[<j)^(z)] 

+ Re [ p ( z ) ] } . 

Ev = (3-v) Im [<}> (z) ] + ( 1+v ) {xlm [<p' (z) ] - yRe [<j> " (z) ] 

+ Im[if> (z) ] } . (12) 

The Goursat-Kolosov equivalent of equations (5) and (7) , with 

* 

the help of (11) , read 

Im[z4>" (z) + \p' (z) ] = 0 

Re[zcj)"(z) + ip'(z)] E A / 

or 

Z(p n (z) + ip' (z) = A (13) 

everywhere. Integrating 

Z<p'(z) - <j)(z) + p(z) = Az + B 

which is equavalent to the pair of equations 

xRe [<p' (z) ] - ylm[4>"(z)] - Re[$(z)J + Re[ijj(z)] = Ax 

yRe t<f> ' (z) ] + xlm[<f> ' (z) ] - Im[<{>(z)] + Im[ty(z)] = Ay 

( 14 ) 
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The constant B, which must be real, can be omitted because its 
retention merely serves to add to the displacement field a 
term which represents a rigid body displacement. Upon combining 
(12) with (14) 

Eu = 2 (1-v) Re [<J> ( z) ] - (l+v)2y Im[(j>''(z)] - (l+v)Ax 

Ev = 4lm[(j) (z) ] - (1+v) 2y Re[<j>"(z)] + (l+v)Ay (15) 


To avoid confusion with the bar symbol used to denote complex 
con jugatic l let 


2<j) (z) 



Z (z) dz 


a- 

Z (z) 


(16) 


where upon 


Eu 

Ev 


'V 

(1-v) Re (Z ( 
a, 


21m [Z (z) ] 


zj] - (l+v)yIm[Z (z) ] 
- (1+v ) y Re [Z (z ) ] + 


- (1+v) Ax 
(1+v) Ay / 


(17) 


emerge as the modified Westergaard field equations for plane 
stress. 

APPLICATIOl fS 

To illustrate use of the modified Westergaard equations 
it is worth, while to treat anew the familiar problem of the 
infinite plate with colinear periodic cracks as shown in Fig. 1. 
The factor k is any real number. 
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Fig. 1. 

Using the Kolosov equations [5] 

a xx + °yy = 2{ ^< z ) + VTZT) = 4Re[<J»'(z)] 

a yy - tf xx + 2ia xy = 2 { z*<f> ” (z) + ip'(z)} (18) 

The boundary conditions can be expressed as follows: For all 


points situated on any crack border 

a yy + ia xy = 2 Re ( <J> ^ ( z ) ] + {z>" (z) + \p'(z)} = 0 ... (19) 

A t | z | = oo 

CT yy (°°) = a xx (oo) = ka ' a xy (oo) = 0 (2°) 
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Due to the symmetry of the loading relative to the x axis (13) 
must be satisfied. With z = z at y = 0, (13) reduces (19) to 

2Re [<p' (z) ] = - A (21) 

for all points on. any crack border. In semi-inverse fashion, 
owing to the periodic and symmetric- nature of the crack spacing, 
the function 2<J> y, (z) can be chosen to have the form 


2<T(z) = 


g (z) 

2,rz. . 2,ira 

sin ‘ sin ( “w 


) 1/2 

>1 


( 22 ) 


where the- denominator of the first term has no real part along 

* 

the crack borders. The function g(z) is presumed to be holc- 
morphic in the region of definition, except possibly at the 
point z = », and must be such that Im[g(x)] = 0 along the crack 
borders. The function 2<j>''(z)' so defined satisfies boundary con- 
dition (21) . 

From boundary condition (20) 


a (°°) - cr (oo) + 2ia (oo) = (l-k)a = 4ylm[<j>" (z) ] + 2A 
jrjf xy 

- 4iy Re [<j) " ( z) ] 

from which 

(1-k) a = (z-z)2<T(z) + 2A, [ z | + «>. .. (23) 

Inasmuch as \2§' {z). must be holomorphic throughout, including 
the point at infinity, it will therefore be continuous at and 
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in the neighborhood of this point, and for |z| arbitrarily 
large 

2 <T(z> - ~ 3 ~ - A 

sin (£L) 

where upon 



which can be identically satisfied by choosing 
g(z) = a sin (~) 

A := j (l-k)a . (25) 

The condition that Im[g(x)] = 0 along the crack borders is 
also seen to be satisfied. 

The stress function which solves this problem is thus 


sin(jp) 

tsin 2 (^) - sin 2 (S.)}' 


For uniaxial uniform tension applied in the y direction, 
k = 0 and A = a/2. The stress function (26) then assumes 
a form equivalent to that given by Sanders [6] . When k = 1, 
A = 0, which corresponds to loading by equal uniform biaxial 
tension. The stress function introduced by Westergaard for 
this problem in reference (3] is therefore a solution only 
for this special loading condition. 
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As another illustration of consequence concerning this 
particular class of plane crack problems, consider the so- 
called crack-tip stress and displacement field equations. 
These can be obtained for opening mode crack surface displace 
ments (mode I) by consideration of the problem of Fig. 1, 
modified to a single centrally located crack of length 2a. 

A stress function which will satisfy the boundary conditions 
along such a cut has the form 


2cf> " (z) = Z (z) = 


g(z) 
2 2 


{z -a } 


1/2 


- A . 


Proceeding as in the previous' example, it will turn out that 
g (z) = crz and a = (l-k)a/2 so that 


Z(z) = 


h 2 -a 2 ) 


£? _ M-v) — 

1/2 u K ' 2 


(27) 


Introducing crack-tip polar coordinates (r,0) through the 


coordinate transformation £ = (z-a) = re 


i0 


Z(?) = 


g (g +a ) _ n_ k , £ 

1/2 ■ u K 2 ' 


{U+a) 2 -a 2 } 


For |?| very small, i.e., |?| <<a 

K. 


z(c) = 


172 " (1 " k) I 

{2tt?} 1/2 2 


(28) 


where 


Kj = 0 {ira} 


1/2 


(29) 
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is the crack-tip stress intensity factor. Substituting (28) into 
(9) , (16) and (17) one obtains for the plane stress crack-tip 
stress and displacement fields the approximations 
K 

V ° I 2J} 3 -/ 2 cos< ^ [1 " sin( ^ sin< i 1)! ' (1 - k)0 


'V 


yy {2iTr} 


K 

— I -y/ 2 cos (|) [1+ sin(-|) sin(|i)] 


K 

a xv * *~ T /? sin cos cos 

Xy {2rrr} 1/2 2 2 2 

K 1/2 

U = / {frF^ cos (|) t(T+^) + sin 2 (|-)]- (l-k)r cos 0 


v = r { f ? )1/2 sin( !» 00s 2 (f> ] 


+ (l-k)r sin 0 .... (30) 


E 


Again only when k = 1, i.e., equal uniform biaxial tensile 
loading, do these equations reduce to the form currently found 
in the literature [7] . 

To further illustrate use of the modified Westergaard 
equations consider the centrally cracked strip (plate) of 
finite width loaded uniaxially in uniform tension, Fig. 2, of 
great interest in fracture toughness testing, and which has 
not been given an exact closed form solution. 



c; 
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Fig. 2 

A widely used approximate solution to this problem was 
first introduced by Irwin [4] , by means of the stress func- 
tion (26) with k = 1, which, as has been shown, is the exact 
solution tc the periodic colinear crack problem in an infinite 
sheet loaded in uniform biaxial tension. To the stress field 
associated with this stress function Irwin adds a uniform 
horizontal compressive stress of magnitude o along the vertical 
edges of the strip which, interestingly, has the effect of 
compensating for the missing A term. This combination satisfies 
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boundary conditions along the vertical edges, leaving a 
horizontal stress of varying magnitude which depends on the 
relative crack size. The crack tip stress intensity factor 
emenating from this stress function is the so-called tangent 
formula 

. ■ . „ 1/2 ' 

K r = o{W tan } ( 3 i) 


Subsequently more accurate truncated series (poly- 
normial) representations for Kj have obtained by Isida [8] 
and Srawley et.al. [9], which show the tangent formula to be 
in varying degree of small error, depending on the crack size. 
Recently a secant formula has been proposed by Fedderson [9] 


K-j. = o' {fra sec 


(— ) } 


1/2 


(32 


/ 

which matches almost identically Isida’ s Kj values, deemed 
to be the most accurate. Having the added virtue of being 
concise and therefore relatively simple to use, Fedderson’ s 
secant formula has now in some quarters replaced the tangent 
formula in fracture toughness testing. 

There will be some practical interest then in obtaining 
the corresponding stress function, that is, one which comes 
acceptably close to solving the problem of Fig. 2 and which 
yields the secant formula for . 

It is convenient to let 


Z(z) 5 ztz) - A . 


!33) 
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Then 


a xx ~ Re [Z* (z) ] - ylm [Z (z) ] - 2A 


a yy = Re [Z* (z) ] + ylm [Z*(z)] 


( 34 ) 


a xy “ y Re ( z )] 


Z (z) 


J 


(Z*(z) - A) dz = Z*(z) - Az 


(35) 


and 


Eu = : (1-v) Re[Z*(z)] - (1+v) yIm[Z*(z)] - ~2Ax 


Ev =■• 2Im[Z* (z) ] - (1+v) yRe [51* (z) ] + 2vAy 


(36) 


A stress function which satisfies the crack border 
boundary condition, partially satisfies the vertical edge 
boundary condition and yields the se:cant formula for Kj has 


the form 


Z (z) = Z* (z) - A = 


1/2 

„ rTra /fra. ' . ,7rz. 

olr-r- csc( r7 -) } sm( fT -) 




'W 


w 


r . 2 /Ttz.- . 2 ,7ra. 

^ sin “ sir ( ir )} 


1 1/2 
1 _ r ^ a ,7Tav T ' 

- 2 a V csc( W _)) 


(37 


For ? << a, where Z = z-a = re 


ie 


z(z) - 


■ fira. #7Tci x 

a{ w~ csc( w _) } 


1/2 


2tt? . ,ira x ,TTa. 

w sxn( ^w~ ] cos( vr ) 


.2, -rra. 2 tt£ . ,7ra. ,ya. 

sin (pT - ) + s i n (vr) cos( r ) 


1/2 
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from which 


* i\j 

z (t) = 


a{?ra sec (—) } 


1/2 


K. 




ITT 


{ 2 *r 5 ) 1/2 


(38) 


Using (38) and 

2A = esc (^) } (39) 

in (34) thru (36) will give the crack tip stresses and dis- 
placements as in equations (30), except that (l-k)a is 
replaced by 2A as given by (39) . 

The crack border condition 


a yy = a X y = °' y = W < * 


is seen to be satisfied by inspection. At z = ~ + iy 


Z (" + iy) 


7T v ■ 


i2A ^ tanh (^-) 


• / n 'a \ — \ 

sxn w ~ 

_COSh(;p^)_ 


1 - 


cosh^j 


-3/2 


which has no real part. Thus a xy (|>y) = 0 for all y. On the 
other hand 


\ 


°xx<!'*> 


rsin(^-) *r r -* 

[ 1+ 1 ' : tanh( ^»] 



- i 


(40) 
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(37) would be an exact solution if the right side of (40) were 
to vanish for all values of y, all other boundary conditions 
having been satisfied. Results of calculation of (40) are 
shown in Fig. 3. For small crack sizes, (ira/W) £ 0.3, the 
right side of (40) gives values values very close to zero 
along the entire vertical edge, having a maximum value of 
about four percent of the applied load at the crack plane 
when (Tra/W) = 0.3. For (ira/W ) > 0.5 the resulting horizontal 
boundary stress exceeds fifteen percent of the applied load 
at the crack plane. The pattern of this boundary stress 
distribution is interesting in that .through Poisson's Ratio 
effects it tends to suppress vertical displacement of points 
situated just above and below the crack plane. 


Owing to the greater relative accuracy of the secant 
formula for K^, one might expect that for small to moderate 
crack sizes, e.g., (ira/W) < 0.5, the stress function (37) 
will yield good estimates for other centrally located quantities 
such as the crack opening displacement, of interest in elastic 
compliance calibrations. For displacement gage points located 
along the plate center line, it follows from (36), after some 
calculation, that 


E 

0 W 


v (o ,y) 


/ 7T 3. / 7T 3. x 

V csc <w- )} 


1/2 


■ 

r - 


c °sh(^) 

— cosh 

■». 

cos (|£) 

L- 




(1+v) 

W 


1 

w — 

2\ 

( 1+ 

sin(£ 1 -) 

1 

sin (^) 

f 

I 


- . -J 

J 


- 1/2 
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Calculation of (41) is compared with experimental data 
obtained from Alum. 7075-T6 center cracked sheets , reported 
in reference [10] , and shown in Fig. 4. The data points 
defining the experimental curve were obtained in the low load 
or elastic range. The predicted crack opening displacement, 
eq . (41) , is a considerable improvement over Irwin's cal- 
culation, and is surprisingly close to the experimental curve 
in the large crack size range where the vertical edge boundary 
condition is poorly approximated. The fact that the predicted 
compliance curve lies entirely below the experimental curve 
appears to be explainable by the particular nature of the dis- 
tribution of the excess of vertical edge boundary stress shown 
in Fig. 3. Imposition of an identical distribution along the 
vertical edges, but reversed in sense, (leaving those edges 
free of traction as they should be) would tend to increase 
somewhat the vertical displacement from that given by (42) for 
all points a little above and below the crack plane and would 
thereby elevate the curve of (41) . 
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